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1. Introduction 

One of the major challenges since the discovery of QCD has been to understand 
hadronic Fock-space wave functions in the high-energy limit. It was realized long 
ago that due to the soft singularity of the splitting function of non-Abelian gauge 
bosons, hadronic wave functions would contain a large number of gluons with a small 
(light-cone) momentum fraction x.^ This expectation was confirmed hy e + p deep- 
inelastic scattering experiments at HERA (see below) and led to the replacement 
of early gluon distributions with little "small- a; glue", such as the Duke-Owens 
distribution with xg{x) w const at — 4 GeV^, by modern PDFs like CTEQ with 
a high density of gluons at low x. 

The growth of the gluon density with decreasing x (or increasing rapidity) im- 
plies that at sufficiently high energy non-linear dynamics of the soft color fields 
would emerge. That is, the increase of the gluon density within a rapidity step from 
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Y to Y + AY is no longer simply proportional to the number of gluons that were 
already present at FpiHor else unitarity is violated eventually. Instead, the scatter- 
ing amplitude is thought to saturate at some dynamically generated virtuality scale 
Qs{Y) which increases with rapidity. This new phenomenon is commonly referred 
to as saturation of gluon distributions at low x and manifests in non-linear, density 
dependent terms in the QCD evolution equations. Furthermore, high gluon densities 
enhance the role of multiparton interactions during the collision process, known as 
higher twists in the context of collinear factorization. At moderate transverse mo- 
mentum multiple scatterings may become a leading effect and should be resummed 
to all orders, which implies a rearrangement of the perturbative series. 

The saturation momentum Qg is boosted further in heavy nuclei due to the fact 
that the density of large- a; "valence" charges per unit transverse area increases in 
proportion to the thickness ~ A-^^^ of a nucleus.!^^ These represent the sources 
for the small- a; soft gluon fields. Thus, non-linear color field dynamics should be 
operative at higher transverse momentum scales than for nucleons. Investigating 
the gluon structure of a hadron or nucleus at low x is of fundamental interest, 
and this can be well done using e-|-p, e-|-A and p-l-A collisions. Understanding the 
dynamics of small-x gluons in nuclei is also a crucial factor in determining the initial 
state of a heavy ion (A-|-A) collision right after impact. This is the aspect which 
we emphasize in this review, with a perspective on the latest results from the LHC 
heavy-ion program. 

Different approaches have been developed to study initial state effects, we shall 
focus on the Color Glass Condensate (CGC) formalism while briefiy discussing 
alternatives in parallel. Knowing the density and distribution of gluons that have 
been released is necessary for a complete dynamical description of the collision 
before the eventual thermalization of the system and the formation of a Quark Gluon 
Plasma (QGP). The practical implications are many. In the soft sector, say particles 
with small transverse momentum pt J; 1 GeV, the initial conditions determine the 
bulk features of multiparticle production in heavy ion collisions such as dN/dy 
or dEr /dy as a. function of collision energy, the collision geometry (distribution of 
produced gluons in the transverse plane) and so on. A very important application of 
the CGC formalism is to provide such initial conditions for the subsequent evolution. 
This input affects significantly the transport coefficients such as shear viscosity 
extracted from hydrodynamic analysis of heavy-ion collisions. Also, it gives rise to 
non-trivial long-range correlations in rapidity and to density "hot spots" which are 
visible in final-state azimuthal correlations. 

In the "hard probes" sector, i.e. for particles with a -perturbatively- large trans- 
verse momentum that do not thermalize but are used as tomographic probes, a 
proper distinction of initial state effects from those due to the QGP ("final state 
effects" ) is vital for a quantitative characterization of the matter produced in heavy 
ion collisions as they may sometimes lead to qualitatively similar phenomena in 
observables of interest. 

Last but not least, another primary goal of studies of the initial state is to 
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provide proof that (local) thermalization of the system actually happens over the 
time scales estimated from hydrodynamical simulations, Ttherm ^ 0.5-^1 fm/c. How 
such ultra-fast isotropization of the system happens is, arguably, one of the most 
fundamental open problems in the field of heavy ions. 

One important lesson learned from experimental data collected in Au+Au and 
Pb+Pb collisions at RHIC and the LHC, respectively, is that bulk particle produc- 
tion in ion-ion collisions is very different from a simple superposition of nucleon- 
nucleon collisions. This is evident, for instance, from the measured charged particle 
multiplicities which exhibit a strong deviation from scaling with the number of in- 
dependent nucleon-nucleon collisions: '^^^ {rj = 0) <C Ncoii '^^^ [r] = 0). One is 
led to the conclusion that strong coherence effects among the constituent nucleons, 
or the relevant degrees of freedom at the sub-nucleon level, must be present during 
the collision process. The fact that coherence effects are essential for the description 
of data is confirmed by the fact that most -if not all- successful phenomenological 
approaches include them in one form or another. In a QCD description, coherence 
effects can be related to the presence of large gluon densities both prior to the 
collision and during the collision process itselo 

2. Color Glass Condensate 
2.1. A Brief review 

The CGC is formulated as an effective theory for the small- a; degrees of freedom 
in high-energy QCD scattering. It relies on a separation of degrees of freedom into 
soft, small- a; gluons -treated as dynamical gauge fields- and high-x, static valence 
degrees of freedom (in terms of light-cone time). The latter act as random sources 
for the small-x gluons. Any physical observable is averaged over all the possible 
configurations of the sources according to a probability distribution W^[p]y. The 
fact that the averaging is performed at the level of the observable refiects the fact 
that at very high energy the hadron is probed in a frozen configuration and, hence, 
interferences between different field configurations are neglected. This is thus an 
essentially classical picture. 

The CGC is endowed with a set on non-linear renormalization group equa- 
tions -the B- JIMWLK equation^sHUL for the probability distribution of the valence 
sources. This ensures that observables do not depend on the particular momentum 
scale at which the separation between soft (dynamical) and valence (static) de- 
grees of freedom is performed. The B- JIMWLK equations account for quantum 
corrections and, at leading logarithmic accuracy, they resum terms ~ aslii(l/a;) 
to all orders. They also include non-linear terms that, in the proper gauge, can be 
interpreted as due to gluon-gluon recombination processes. The mere presence of 
non-linear terms immediately implies that a dynamical transverse momentum scale 
emerges, the saturation scale Qs, such that gluon modes with transverse momentum 
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kt ^ Qs{x) are in the saturation regime. 

Saturation of gluon densities is equivalent to the presence of strong color fields, 
parametrically of the order of the inverse of the coupling, A{x) ^ 1/g. Thus, al- 
though the insertion of new sources in diagrams for particle production appears 
to be suppressed by powers of the coupling constant, this is compensated by the 
strength of the color fields, i.e. terms of the order g A{x) ~ 0{1) must be resummed 
to all orders. This implies a rearrangement of the perturbative series in the high- 
density regime or equivalently, to an all-order resummation of multiple scatterings. 
The propagation of an energetic projectile, e.g. a quark, moving along the (light- 
cone) direction, at some transverse position z± in the background of the strong field 
of a nucleus, denoted A, is conveniently described in terms of Wilson lines by (we 
are using the A~ = gauge): 



Wilson lines resum multiple scatterings in the eikonal approximation (neglecting 
the recoil of the projectile) and they are the relevant degrees of freedom in the 
high energy regime, rather than free quarks or gluons. Actually, the B-JIMWLK 
equations can be recast as an infinite hierarchy of coupled, non-linear differential 
equations for the rapidity {Y = ln(l/a;)) evolution of n-point correlators of Wil- 
son lines, i.e. for correlators of the type {U{x±i) . . .U{x±n))Y ■ The B-JIMWLK 
equations are difficult to solve numerically and to implement in phenomenological 
applications. In turn, the BK equation,'i2l!ill which corresponds to the large- A^c limit 
of the full B-JIMWLK hierarchy and provides a single, closed evolution equation 
for the 2-point function, {U{x±i)U^ {x±2))y , is much easier to solve in practice and 
is widely used in phenomenological applications. As we shall discuss in detail below, 
the precise prescription for the resummation of multiple interactions varies depend- 
ing on the observable of interest or the colliding system, and does not always result 
in easy-to-use factorizable expressions. 

A different approach to particle production at high transverse momentum as- 
sumes coUinear factorization of heavy ion collisions as a working hypothesis. In 
this framework all nuclear effects -and in particular nuclear shadowing at small- cc- 
are absorbed into modified nuclear parton distribution functions (nPDF's) fitted 
to experimental data; the scale dependence is described by DGLAP evolution. A 
variety of intermediate approaches can be found in the literature that start from 
the coUinear factorization ansatz and further extend it to account for multiple scat- 
terings through the calculation of higher twists, in the coherent case, or through an 
all order resummation of incoherent, Glauber like scatterings. The result of these 
calculations is sometimes recast in the form of modified nPDF's which include an 
energy-dependent intrinsic transverse momentum and momentum broadening. Fur- 
ther, cold nuclear matter energy loss is sometimes included in these formalisms. 
These essentially perturbative tools are complemented with non-perturbativc mod- 
els for soft particle production in Monte Carlo event generators such as HIJINGf^or 
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HYDJET.^^ There, phase-space arguments motivate the implementation of energy- 
dependent transverse momentum cut-offs (analogous, in essence, to the CGC sat- 
uration scale) in order to regulate independent particle production from different 
sources in the hard sector. A detailed compilation of such different approaches can 
be found inl^^ 



2.2. Theory Status 

Despite the intrinsic difficulty of performing perturbative calculations in the back- 
ground of a strong color field, steady progress in the degree of accuracy of the CGC 
formalism has been achieved over the last years. Briefiy, it can be said that the 
CGC is now entering the next-to-leading (NLO) order era. 

Thus, the kernel of the BK equation is now known to NLO accuracy in the resum- 
mation paiametei^ as \n{l/x) or, through an all-orders resummation of a partial 
subset of NLO corrections, to running coupling accuracy.^^^^ The BK equation 
including running coupling corrections (henceforth referred to as rcBK equation) 
reads 



|r^^= / rf'li i^™"(z:,ri,r2)[A/>(ri,a;)+AAf(r2,a;)-AA;^(r,x) 
dln(a;o/a;) J 

-NF{rux)NF{r2,x)] , (2) 

where J\fF{r = \z_j^ — z_2\,x) — l — ^{tr [U{z^i)U^ {z^2)])Y=ini/x is the (imaginary part 
of) dipole scattering amplitude in the fundamental representation, and = r — r^. 
xq is some initial scale and K^^" is the evolution kernel including running coupling 
corrections. Different prescriptions have been proposed in the literature for i^f"". 
As shown in ref.l^ Balitsky's prescription minimizes the role of higher conformal 
corrections: 



i^™"(r,ri,r2) 



27r2 



2\ 



(3) 



Analogous effort has been dedicated to improve the degree of accuracy of the ex- 
pressions for particle production processes. Thus, a full NLO generalization of the 
dipole factorization formula for the calculation of structure functions in deep in- 
elastic scattering (DIS) at small-x is now available .l^l^ It should be noted, though, 
that the AAMQS fits reported below were performed using the LO dipole formal- 
ism. Also, corrections (sub-leading in density) to the McLerran-Venugopalan action 
that provide a theoretical justification for the AAMQS parametrization of the initial 
conditions of the evolution (see below) have been recently calculated.'^ 

In the case of proton-nucleus collisions there are two distinct but related ap- 
proaches to single hadron production in the CGC framework: fci-factorization and 
the hybrid formalism. Their applicability depends on the kinematic region in which 
the projectile proton is probed: small-x or high-x respectively. In both cases explic- 
itly factorized expressions in terms of only the 2-point function or, equivalently, the 
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unintegrated gluon distribution, are well established. Actually, running coupling, 
and full NLO corrections to the /ci-factorization and hybrid formalisms have been 
calculated iiP ancP^ respectively. However only the inelastic corrections -a sub- 
set of the full NLO corrections- to the hybrid formalism calculated irP^ have been 
implemented in phenomenological works. These first studies on the impact of 
the inelastic correction indicate that they are sizabl^^ and, in some cases, they 
overwhelm the LO contribution in certain kinematic regions (high pt), which calls 
for a phenomenological study of the full NLO corrections. 

The knowledge of the 2-point function provided by the BK equation docs not 
suffice to describe more exclusive observables. Rather, knowledge of higher point 
functions, and hence of solutions of the B-JIMWLK equations are needed. An ex- 
tensive study of the different correlators involved in the calculation of two-particle 
production processes (di-jet or photon-jet) in p-l-A collisions and DIS was presented 
injm While the techniques for solving the B-JIMWLK equations numerically have 
been known for some time, their application remains somewhat difficult. Significant 
progress has been made recently in devising approximate analytical solutions for 
arbitrary n-point functions based on the Gaussian truncation of the B-JIMWLK 
hierarchy,'^^t23i ^Ymg extending the semi-classical MV model for the initial condi- 
tion to evolved systems. This new technique allows to express arbitrary n-point 
functions in term of the BK-evolved 2-point function thus opening new avenues for 
a systematic phenomenological exploration of multigluon correlations in available 
data. A first application of the Gaussian approximation to the study of di-hadron 
correlations at RHIC can be found in '^'^ where the role of double parton interactions 
-important for the analysis of di-hadron correlations data-in the CGC has also been 
investigated. 

In the more intricate case of nucleus-nucleus collisions all leading terms in the 
calculation of inclusive particle production can be resummed by solving the classi- 
cal Yang-Mills equations of motion (GYM), with the valence degrees of freedom of 
the two colliding nuclei acting as external non-dynamical sources. However, supple- 
menting the classical methods with information on the quantum evolution of the 
wave function of the colliding nuclei, necessary to extrapolate from RHIC to LHG 
energies, is presently a difficult task since, again, it requires the numerical resolu- 
tion of the full B-JIMWLK equations. This difficulty is bypassed in the IP-Glasma 
model, '^^ where B-JIMWLK evolution is replaced by the empiric parametrization of 
the saturation scale obtained through the IP-Sat model fit to HERA data. Neverthe- 
less, a proof of factorization of inclusive gluon production and multigluon correla- 
tors including quantum corrections has been provided in.lSfill^Z! Those results provide 
a dynamical framework for the Glasma flux tube picture and provide a rigorous 
theoretical basis for the study of long range in rapidity correlations in the CGG 
framework, as we discuss below. 

Finally, intense effort has been devoted to the problem of thermalization in heavy 
ion collisions in the recent years. It has been recently demonstrated in ref.'^* that 
the rcsummation of the leading secular terms in the evolution of a scalar 0^ theory 
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coupled to a strong external source makes the pressure tensor of the system relax 
to its equilibrium value, a prerequisite for the onset of hydrodynamical flow. An 
equivalent result for QCD would provide a proof of thermalization of the system 
and may allow to determine the equilibration time as a function of the natural 
scales of the problem (collision energy, energy density etc) . For a recent review on 
this topic see, e.g.!^ 

Finally, it is worth remarking that the substantial theoretical progress reported 
above has not been yet fully exploited in phenomenological applications, which 
leaves a large margin for improvement in the analysis of the upcoming LHC data. 

3. e+p to e+A collisions 

Despite the fact that CGC effects are expected to be enhanced in nuclei versus 
protons, which is due to the larger valence charge densities per unit transverse area 
in nuclei, so far the most exhaustive searches for the gluon saturation phenomenon 
have been performed using data on proton reactions. This is mainly due to the large 
body of high quality DIS data on protons at small-x. Moreover, most -if not all- 
phenomenological approaches for heavy ion collisions borrow empiric information 
from the analysis of e+p data, either to parametrize the Bjorken-a; or energy de- 
pendence of the saturation scale or to also constrain more exclusive features of the 
unintegrated gluon distribution of the proton. 

A first analysis of DIS data in terms of saturation physics was provided by the 
GBW model.^ There, a good description of data was obtained with a saturation 
scale QsN(a;) = (^) ' GeV^ with xq ~ 3 ■ 10^'*. Such empiric parametrization 
found a theoretical justification in the analysis of the collinearly improved Leading 
and Next to Leading BFKL evolution in the presence of saturation boundaries 
performed in.lSl Other analytic models that explicitly include saturation or unitarity 
constraints and which also incorporate features of BFKL dynamics, like the IIM 
niodel,'^ or are based on eikonalization of leading twist scattering in the collinear 
factorization formalism like the IP-Sat model,'^^ also provided a good description of 
the dat^EI and are frequently used for phenomenology. 

An important step forward towards a theory driven analysis of the data was 
brought about by the AAMQS global fit^^sillfii inclusive structure functions mea- 
sured in e+p. The AAMQS approach relies on the LO dipole model plus the rcBK 
evolution equation to describe the small-x dependence of the dipole scattering am- 
plitude, Af{r, x). The free parameters in the AAMQS fits, aside from a global nor- 
malization constant and an infrared parameter related to the argument of the run- 
ning coupling (se^^ for full details) , correspond to the free parameters in the initial 

'^For a recent fit of the IP-sat model to the final data from HERA, including some processes which 
probe the impact parameter dependence of the dipole scattering amplitude, we refer to ref.E^l 
Such modern models do a good job of describing the growth of the gluon density towards low x 
mentioned already in the introduction. 
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conditions for the evolution, as follows: 
7V(r, a; = xo = 0.01) = 1 - exp 



(^2 n2 



'5s0, proton) i / 1 

-r In h e 



Ar 



(4) 



Here, A = 0.241 GeV, Qso, proton is the saturation scale at the initial xq and 7 is a 
dimensionless parameter that controls the steepness of the unintegrated gluon dis- 
tribution for momenta above the initial saturation scale fct > Qso- With this set up, 
the AAMQS fits provide a remarkably good description of data, see Fig. ([T]) right. 
Although the AAMQS fits clearly favor values 7 > 1, they do not uniquely determine 
its optimal value (and neither does the analysis of forward RHIC data performed 
in ref!S3). Rather, different pairs of (Q^q pj.Qto,j, 7)-values provide comparably good 
values of /d.o.f ~ 1, see Table[Tl The "degeneracy" is due to correlations among 
parameters and also because HERA data is too inclusive to constrain exclusive fea- 
tures of the proton UGD. A recent comparative study of the stability of standard 
DGLAP and rcBK fits to small- a; HERA data under changing boundary conditions 
performed by systematically excluding subsets of data from the fits showed that 
rcBK fits do provide more robust fits to data.f^Sl This is an additional indication for 
the relevance of non-linear small-x dynamics in the analysis of HERA data. 



Set 


OsO.proton (GcV^) 


7 


MV 


0.2 


1 


h 


0.168 


1.119 


h' 


0.157 


1.101 



The main obstacle to extending the above studies to the case of nuclear DIS 
is the paucity of experimental data at small-x. This lack of empiric information to 
constrain the non-perturbative parameters in either the collinear factorization or 
CGC formalisms (mostly related to the initial conditions for the evolution) is one 
of the main sources of uncertainty in nPDF parametrizations (see Fig. [T] right) and 
in CGC predictions for the LHC. It should be noted that for nuclear targets one 
additional variable is needed to properly characterize the wave function: the impact 
parameter dependence. This difficulty is bypassed in standard nPDF minimum bias 
fits. A first 6-dependent nPDF set has been proposed recently in.**^ In CGC ap- 
proaches one normally relates the b-dependence of the nuclear saturation scale to 
the proton saturation scale and the local nuclear density: (3^^(6) =rA(&) Q^s,proton- 
This can be done at different levels of refinement: assuming a homogeneous nu- 
cleus, a nucleus with a Woods-Saxon density distribution without fluctuations, or, 
finally, also accounting for geometry fluctuations as in the KLN-MCPIMIand rcBK- 
MCP^I^ Monte Carlo's. Such ansatz, when combined with the AAMQS parameters 
7 > 1 (see Table [Ij may pose problems as it violates the additivity of nuclear dipole 
amplitudes on the number of nucleons at small dipole sizes. In order to fix this 
problem, alternative dependencies of (5^^(&) on TA{h) were explored in ref.l^Data 
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from tht rPb run at the LHC will provide important input for sorting out such 
uncertainties and to get a better determination of nuclear wave functions. 
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- Q^=45 GeV^ ■ 




Fig. 1. Left: rcBK fits to HERA data on the reduced cross section in e+p scattering 
from .EH Right; Different nPDF parametrizations for the nuclear gluon distributiorj^— at 
GeV2. 



at small- a; 
Q2 = 1.69 



4. p+A collisions 

Proton nucleus collisions provide a better handle on initial state and gluon satura- 
tion effects than those of heavy ions. This is mainly due to the absence of final state 
effects induced by the presence of a QGP. There are two distinct but related ap- 
proaches to hadron production in high energy asymmetric (such as proton-nucleus 
or very forward proton-proton and nucleus-nucleus) collisions. In such collisions, 
particle production processes in the central rapidity region probe the wave function 
of both projectile and target at small values of x and can be described in terms of 
the fcf-factorization formalism where the differential cross-section for inclusive gluon 
production in A-|-B collisions is given 



dy (Ppt (PR 



as{Q) 



\Pt 



h\ Apt 
-,xi;b)if{- 



-,X2;R~ b) , 



(5) 



2 '---'-^ 2 

where y, pt and R are the rapidity, transverse momentum and transverse position 
of the produced gluon, Xi(^2) — {pt/ \/snn) exp(±y) and Cf = {N^ — l)/2Nc. The 
relevant UGD's for the fc^-factorization formula above are given by 

Cf 



ip{k^ X, R) 



-zkr 



V^ACi(r,a;,R) 



(6) 



a,(fc) (27r)3 

where Ma is the dipole amplitude in the adjoint representation which, in the large- 
A^c-limit, is related to the dipole amplitude in the fundamental representation by 
AAA(r,a;,R) = 27VF(r, x, R) - A/'|(r, a;, R) 
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However, at more forward rapidities the proton is probed at larger values of 
X while the target nucleus is deeper in the small- a; regime. In that context kt- 
factorization fails to grasp the dominant contribution to the scattering process. 
Rather, the hybrid formalism proposed in^"* is more appropriate. In the hybrid for- 
malism the large- a; degrees of freedom of the proton are described in terms of usual 
parton distribution functions (pdf 's) of collinear factorization with a scale depen- 
dence given by DGLAP evolution equations, while the small- a; glue of the nucleus is 
still described in terms of its UGD. At leading order, inclusive particle production 
is given by: 



72 



drjd^k {2ny 



. 1 



Xlfq/p{xi,Q^)NF ('^2,^) Dh/q{z,Q^) 

xifg/pixi,Q^)NA{x2,^'j Dh/g{z,Q^) , (7) 



where 



Np^A) (fc, R) = Jdh e-^'^-'- [1 - Afp^A) {r, x, R)] . (8) 

Next-to-leading corrections to Eq. ([7]) have been calculated in.!^!^ Clearly, the hy- 
brid formalism is restricted to production of particles with high enough transverse 
momentum, as implied by its use of collinear factorized pdf's and fragmentation 
functions. Unfortunately, no smooth matching between the fct-factorization and hy- 
brid formalisms is known to date. Also, their corresponding limits of applicability 
-equivalently the precise value of x at which one should switch from one to the 
other- have only been estimated on an empirical basis. 



4.1. Transverse momentum integrated yields 

The pT-integrated multiplicity of charged particles is the most basic observable 
in particle collisions. For a variety of reasons though, it is also hard to compute 
with precision, unfortunately. To name a few: in practice, many computations of 
dN/dy rely on fcT-factorization, following the early papers by Kharzeev, Levin and 
NardP^'^ who showed that it works surprising well to describe the centrality de- 
pendence of hadron production even in nucleus-nucleus collisions at RHIC energy, 
although it should not.^^ In another approach, the multiplicity is computed nu- 
merically from an exact solution of the classical Yang-Mills equations, using a MV- 
model action with an energy dependent saturation momentum.'^fiJ This accounts for 
multiple scattering both on the projectile and target side but assumes that the ef- 
fective action far from the valence sources remains a local Gaussian. We expect, 
however, that these calculation will be improved in the foreseeable future to employ 
an effective action that solves JIMWLK. We review hadron multiplicities in A+A 
collisions in section [5731 

The theoretical foundation for fcy-factorization is somewhat more solid in case of 
p+A collisions since dN/dy is dominated by Qs.a ^ Pt ^ Qs.p.'^ Nevertheless some 
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hurdles remain, not the least of which is, of course, hadronization. All CGC based 
phenomenological computations of dN/dy assume that the p^-integrated hadron 
yield is proportional to that of gluons. The number of hadrons per produced gluon 
could depend on energy, though.^^ Furthermore, it is clear that the distribution of 
hadrons in rapidity need not be identical to that of gluons, especially when their 
density per unit of rapidity is high. More to the point: pseudo-rapidity distributions 
dN/drj of unidentified charged hadrons involve a transformation from y ^ rj which 
depends on the mass and pT of a hadron; {pt) increases with energy, with the mass 
numbers of projectile and target, and is affected by soft final-state interactions (at 
least in case of A+ A collisions) . 

It is clear that computing dN/drj distributions from first principles is a 
formidable problem. That said, with some level of optimism one may hope that 
the main dependence of dN/ dr] on energy and system size is through the saturation 
scale Qs- Indeed, one very important outcome of the p-|-A and A-l-A programs at 
RHIC and LHC, in terms of "soft" hadron multiplicities, is that this assumption 
works to a large degree. Thus, the main features of inclusive particle production can 
be understood within weakly-coupled, semi-hard (albeit non-linear) QCD although 
genuine non-perturbative effects such as hadronization physics or soft final-state 
interactions are obviously not entirely irrelevant. 
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Fig. 2. Left: charged particle rapidity distribution for minimum-bias p4-Pb collisions at 5 TeV 
from the KLN gluon saturation model.ES Center: rcBK UGD.^* Right: IP-Sat UGD.iU ALICE 
data from ref.EH 



We now proceed to show some of the most recent computations, within various 
implementations of the CGC formalism, of multiplicities in p-|-Pb collisions at LHC 
energies. We stress that all of these had been made public before the p-|-Pb pilot run 
at the LHC occured in faU 20lfl Fig.[2]shows the prediction of the KLN modeP^lfor 
5 TeV. The original prediction depicted by open squares employed the exact same 
y — !■ 77 transformation, specifically the same rapidity independent ratio of hadron 



Albeit, in most cases the initial predictions assumed a collision energy of y/s = 4.4 TeV which 
was later updated to 5 TeV; also, initial predictions did not account for the rapidity shift by 
Ay = ^log(82/208) ~ —0.465 of the experimental detectors towards the Pb beam and were 
updated later. 
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mass to transverse momentum m/pT — 0.24, as for p+p collisfons. At 77 = the 
prediction essentially coincides with the data, in fact much better than its true level 
of accuracy. Away from midrapidity it exhibits a slightly steeper dependence on rj 
than the data. At such level of detail, however, one can not discard hadronization 
effects entirely; this is illustrated in the figures by the "updated curve" (stars) which 
introduces a slight rapidity dependence of m/{pT), with a sign that follows from 
the rapidity dependence of (pr) ^ Qs.Ph- 

A similar prediction from fcr-factorization with rcBK UGDs is shown in fig. [2] 
(center panel). Once again, dN/drj at rj — matches the data at a level of a few 
percent but the j] dependence predicted with the j/ — > r/ transformation from p+p is 
somewhat too steep. A curve corresponding to rapidity dependent hadron transverse 
momentum px = P + AP{i]) (open triangles) illustrates the sensitivity to detailed 
properties of the hadrons in the final state. 

Yet another prediction for dN/dr] in p+Pb collisions employs fc^-factorization 
with UGDs from the IP-sat model.'^ An updated result for y/s = 5 TeV, and in- 
cluding the above-mentioned rapidity shift, is shown again in fig. [5] on the right. 
These authors also provide uncertainty estimates shown by the bands. This pre- 
diction is also close to the data although it overpredicts the multiplicity at 77 < 1 
somewhat. In part, this may be due to neglecting target thickness fluctuations which 
effectively suppress the multiplicity in the dense regime at pT < Qs (see fig. 7 and 
related discussion in ref.l^Sl). 

Last but not least, Rezaeian provided another set of predictions for ^/s = 4.4 TeV 
in ref."^ which is slightly higher than the above (again, without accounting for target 
thickness fluctuations). In all, it is quite remarkable that all of these implementations 
which differ in their hadronization prescription or in the treatment of the nuclear 
geometry and of Glauber fluctuations lead to similar predictions at 77 = (deviations 
from the central value are ±15% or less); they all involved an extrapolation from 
d+Au collisions at RHIC by a good order of magnitude (factor of 25) in center- 
of-mass energy. This is consistent with expectations that the energy dependence 
of particle production is determined mainly by the growth of the gluon saturation 
scale rather than by genuinely non-perturbative soft physics. 

4.2. Single inclusive hadron yields and nuclear modification 
factors 

Nuclear effects on single particle production are typically evaluated in terms of 
ratios of particle yields called nuclear modification factors: 

, dNP^^'^^/dydW .q. 

N,oU dNPP^hX/dyd^p^ ' 

where Ncoii is the number of binary nucleon-nucleon collisions in a p-|-A collision at a 
given energy. If high-energy nuclear reactions were a mere incoherent superposition 
of nucleon-nucleon collisions then the observed RpA would equal unity. On the 
contrary, RHIC measurements in d+Au collisiond^l^ show strong nuclear effects, 
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with mainly two opposite regimes: at mid rapidities the nuclear modification factors 
exhibit an enhancement in particle production at intermediate momenta \p±\ ~ 2-^4 
GeV. In turn, a suppression at smaller momenta is observed. However, at forward 
rapidities the Cronin enhancement disappears, turning into an almost homogeneous 
suppression. 

Single-hadron production data at mid-rapidity have been successfully analyzed 
through different formalisms and techniques, such as leading-twist perturbation 
theory, Glauber-eikonal multiple scatterings, or the Color Glass Condensate. These 
different approaches offer a comparably good description of the mid-rapidity RHIC 
data, so it is difficult to extract any clean conclusion about the physical origin of 
the Cronin enhancement. This is probably due to the kinematic region probed in 
these measurements. For a hadron momentum of |p_l| ~ 2 GeV, one is sensitive to 
X ~ 0.01^0.1. In this region different physical mechanisms concur so that neither of 
the physical assumptions underlying the various approaches above is fully satisfied. 
Indeed, the coherence length, estimated to be Ic ~ l/{2mNx) '--^ 1 10 fm in the 
target rest frame is, on average, smaller than the nuclear diameter; hence the fully 
coherent treatment of the collision implicit in the CGC formalism is not completely 
justified and large- a; corrections are expected to be relevant. Moreover, both the 
coherent or incoherent treatments of the collision need to invoke the presence of 
an intrinsic scale, presumably of non-perturbative origin, of the order of 1 GeV in 
order to reproduce the data, whose dynamical origin is not evident at all. 

The situation at forward rapidities is better understood. Denoting p±^ and y the 
transverse momentum and rapidity of the final state particles, the partons that can 
contribute to the cross section have a fraction of longitudinal momentum bounded 
from below by Xp (for partons from the proton wave function) and xa (for partons 
from the nucleus wave function) given by 

Xp — xf , XA = Xp e~^^ . (10) 

Here we introduced the Feynman variable xp — \p i]e.y / y/ s n n with ^Jsnn the colli- 
sion energy per nucleon. Therefore with ^/snn S> \p±\ and forward rapidities 2/>0, 
the process features Xp < 1 and x^i <C 1, meaning that forward hadron produc- 
tion is sensitive only to high-momentum partons from the proton (well understood 
in QCD), while only small- momentum partons from the nucleus contribute. This is 
the ideal configuration for probing parton saturation effects, and indeed, saturation- 
based approaches were the only ones to correctly predict the suppression of forward 
particle productior(S2ll22l (as well as the azimuthal decorrelation of forward hadron 
pairs,'^ described in the next sub-section). 

After the forward-rapidity RHIC data confirmed the CGC expectations, it has 
been argued that the observed suppression of particle production at forward ra- 
pidities is not an effect associated to the small values of xa probed in the nuclear 
wave function, but rather to energy-momentum conservation corrections relevant 
for Xi? — )■ l.ESlEa] Such corrections are not present in the CGC which relies on 
the eikonal approximation (this may explain why a iiT-factor is needed to describe 
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the suppression of very forward pions). The energy degradation of the projectile 
parton, neglected in the CGC, through either elastic scattering or induced gluon 
brehmsstralung would be larger for a nuclear than for a proton target, on account of 
the stronger color fields of the former, resulting in the relative suppression observed 
in data. LHC data are crucial to identify the correct suppression mechanism since 
LHC is able to probe small xa while at the same time keeping xp far from the 
kinematic limit so that energy-loss corrections should be negligible. 

Before presenting theoretical predictions and comparisons to LHC data, let us 
discuss the kinematics in more detail. Contrary to popular belief, single-inclusive 
hadron production at mid-rapidity at the LHC does not probe the same values of 
XA in the nucleus as forward hadron proudction at RHIC. The minimal values of 
Xa are indeed similar. However, the mean value is roughly given by xa/{z), where 
(z) is the average momentum fraction carried by a final state hadron relative to its 
mother parton, due to the fragmentation process. At forward rapidity and RHIC 
energy, (z) is close to 1 but at mid-rapidity at the LHC it is around 0.3, meaning 
that larger values of x are probed there. In addition, for p^'^'^^'o" ^ one has 
^parton ^ Q^/(^z), thcsc two cffccts combiuc so that, for similar hadron transverse 
momenta, the suppression of RpA should be less important at mid-rapidity at the 
LHC as compared to what has been measured at forward rapidities at RHIC. 

For 1 < < 3 GeV, the transverse momentum range for which those data 
overlap, this is indeed what CGC calculations predict. This can be seen in the left 
panel of Fig. [S] the CGC predictions are around RpPb ~ 0.6-;- 0.9, and within errors 
they agree with recent ALICE data. At higher pr the CGC predictions are also 
confirmed but we emphasize that this would not have been the case without some 
recent development: the predictions for instance are well below the ALICE 
data. The main reason is that, in this latter study and in other earlier works, the 
impact-parameter dependence of the nuclear gluon density, and of the saturation 
scale, has not been described properly: minimum-bias predictions were obtained 
using an impact-parameter averaged saturation scale. Later CGC studies such a^^H 
have shown that it is very important to account for the proper nuclear geometry as 
well as for its attendant fluctuations. 

Other theory comparisons with the ALICE mid-rapidity data show that, other 
than parton saturation effects, there are no cold matter effects. In particular there 
is no Cronin enhancement at low px and no suppression at high pT- For pT > 4 
GeV, RpPb is compatible with unity and well described by the leading-twist nuclear 
PDF approximation. Further data at backward rapidities, probing larger values of 
Xa, will be needed to confirm the anti-showing peak predicted in this approach. 
Conversely, additional data for forward rapidities, probing smaller values of xa, 
will be needed to confirm the suppression pattern predicted in the CGC approach, 
as shown in the right panel of Fig. [3l 
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Fig. 3. Left: ALICE preliminary datsPH for the charged hadron nuclear modification factor at 
m id-rap idity as a function of pT- The theoretical results correspond to the CGC calculations 
ofpSmiEU the nuclear PDF approach EPS09,^* the cold-nuclear matter predictions of, -6^^ and the 
HIJING Monte Carlot^ll Right: LO EPS09^ and rcBK-MO^Si forward-rapidity predictions. 



4.3. Double inclusive particle production and azimuthal 
correlations 

We conclude this section on p-l-A collisions with a discussion of forward di-hadron 
correlations measured at RHIC (ridge-like correlations of rapidity separated di- 
hadrons measured at the LHC are discussed later in section [5^ . In the case of 
double-inclusive hadron production pA^hih2X, denoting pi±, P2± and j/i, y2 the 
transverse momenta and rapidities of the final-state particles, the Feynman variables 
are Xi = \pi±\e^^ / ^/snn and Xp and xa read 

Xp = xi + X2 , XA ^ xi e^^^^ + X2 e^^^^ . (11) 

We shall only consider here the production of two forward particles, since this is the 
only case which is sensitive to values of x as small as in the single- inclusive case: Xp < 
1 and XA^i- The central- forward measurement does not probe such kinematics: 
moving one particle forward increases significantly the value of Xp compared to 
the central-central case (for which Xp = xa = \p±\/ ^/snn), but decreases xa only 
marginally. In addition, we will focus on the A0 dependence of the double-inclusive 
hadron spectrum, where Acf) is the difference between the azimuthal angles of the 
measured particles hi and /i2- 
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The kineraatic range for forward particle detection at RHIC is such that Xp^OA 
and Xyi~10~^. Therefore the dominant partonic subprocess is initiated by valence 
quarks in the proton and, at lowest order in as, the pA^ hih2X double- inclusive 
cross-section is obtained from the qA — >■ qgX cross-section, the valence quark density 
in the proton fq/p, and the appropriate hadron fragmentation functions -D/i/g and 

Dh/g- 

1 rl 



dNP^^'''''-''{P,p,,P2)^ I dzi dz2 I ^ dx fq/p{x,fi^) 

^^,A^,gX Lp Pl PA DhJq{z^,fi^)Dh,/g{z2,fl' 



Zl Z2 



^^.A^ggX Lp El EA Dh,/g{zi,f,')Dh,/q{z2,fl') 
\ Z2 Z\ J 



(12) 



We have assumed independent fragmentation of the two final-state hadrons, there- 
fore Eq. (fT2|) cannot be used to compute the near-side peak around A$ = 0. Doing 
so would require the use of poorly-known di-hadron fragmentation functions, rather 
we will focus on the away-side peak around A$ = tt where saturation effects are 
important. Note also that if Xp would be smaller (this will be the case at the LHC), 
the gluon initiated processes gA qqX and gA ggX should also be included in 
Eq. (ini). 

On the nucleus side, the gluon longitudinal momentum fraction varies between 
xa and e^^^^ -|-e^^^^. Therefore with large enough rapidities, only the small— a; part 
of the nuclear wave function is relevant when calculating the qA — > qgX cross section 
(there is no contamination from large x components). When probing the saturation 
regime, one expects dN'^^^'^^'^ to be a non-linear function of the nuclear gluon 
distribution, i.e. that this cross section cannot be factorized further: dN'^^'^'^^^ ^ 
fg/A ^ dN'^B^'^s^ . Using the CGC approach to describe the small— cc part of the 
nucleus wave function, the qA -> qgX cross section was calculated in,'^'^2Hli] and 
indeed it was found that, due to the fact that small-x gluons in the nuclear wave 
function behave coherently and not individually, the nucleus cannot be described 
by only a single-gluon distribution. 

The qA^qgX cross section is instead expressed in terms of correlators of Wilson 
lines (which account for multiple scatterings), with up to a six-point correlator 
averaged over the CGC wave function. At the moment, it is practically difficult 
to evaluate the six-point function. In the large- A'^c limit, it factorizes into a dipole 
scattering amplitude times a trace of four Wilson lines (or quadrupole), and the 
latter can in principle be obtained by solving an evolution equation written down 
in.f^ However, this implies a significant amount of numerical work, and also requires 
to introduce an unknown initial condition. In practice, further approximations have 
been used, and this is where different phenomenological studies of forward di-hadron 
correlations in p-l-A collisions differ: 



• in Ref.,'^the quadrupole is simply disregarded, the cross-section is obtained 
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solely from the dipole amplitude, using the so-called /cr-factorized formula 
recovered in the limit Qs/\p±i.2\ 1- Saturation effects are nevertheless 
included in the UDGs, and it is done using the GBW parametrization!^ 

• in Ref.,l^the quadrupole is evaluated using the so-called Gaussian approx- 
imation of B-JIMWLK evolution, however only the elastic contribution is 
kept. The non- linear evolution is obtained using the rcBK equation. 

• in Ref.,'^the complete Gaussian expression of the quadrupole is used, how- 
ever only the so-called correlation limit Qs ^ \p±i +P±2\ ^ |Pii,2| is con- 
sidered. Saturation effects are included using the GBW parametrization. 
The gluon-initiated processes calculated irPslare included for the first time. 

• in Ref.|^ the complete Gaussian expression of the quadrupole is used, 
and the non-linear evolution is obtained using the rcBK equation. Gluon- 
initiated processes are not included. 

We now come to the comparison with data. Nuclear effects on di-hadron cor- 
relations are typically evaluated in terms of the coincidence probability to, given 
a trigger particle in a certain momentum range, produce an associated particle 
in another momentum range. The coincidence probability is given by CP{A(f>) = 

Npair{A(l))/Ntrig with 




-1 1 2 3 4 5 



ACP 

Fig. 4. STAR preliminary data^^ for the Coincident Probability between pairs of hadrons as 
a function of the relative azimuthal angle in d-l-Au collisions at RHIC. The theoretical results 
correspond to two CGC-based calculation^Sll an<iS31 and a higher-twist one.ESl Figure fronJSSl 



The STAR data for the coincidence probability obtained with two neutral pions 
are displayed in Fig. 2] for central d-f Au collisions. The nuclear modification of the 
di-pion azimuthal correlation is quite impressive, considering that the prominent 
away-side peak seen in p+p collisions is absent in central d-|-Au collisions, in agree- 
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ment with the behavior first predicted inP^ The data are compared with the CGC 
calculations oP^ andj^ as well as with a non-CGC approach discussed below. As 
mentioned before, the complete near-side peak is not computed, as Eq. ([T2|) does 
not apply around = 0. Note also that, since uncorrelated background has not 
been extracted from the data, the overall normalization of the theory points has 
been adjusted by adding a constant shift. 

The physics of the disappearance of the away-side peak is the following. The 
two measured hadrons predominantly come from a quark and a gluon, which were 
back-to-back while part of the initiating valence quark wave function. During the 
interaction, if they are put on shell by a single parton from the target carrying zero 
transverse momentum, as is the case when non- linear effects are not important, then 
the hadrons are emitted back-to-back (up to a possible transverse momentum broad- 
ening during the fragmentation process). By contrast, in the saturation regime, the 
quark and antiquark receive a coherent transverse momentum kick whose magni- 
tude is of order Qs, which depletes the correlation function around A(f> = n, for 
hadron momenta not much higher than Qs- 

The sizable width of the away-side peak in Fig. 2] cannot be described within 
the leading-twist collinear factorization framework, which completely neglects non- 
linear effects. While it may not be obvious from single-inclusive measurements, this 
indicates that power corrections are important when ~ 2 GeV. At such low 
transverse momenta, collinear factorization does not provide a good description 
of particle production in p-f A collisions. In contrast, CGC calculations which do 
incorporate non-linear effects by resuming power corrections can reproduce the sup- 
pression phenomenon. Ref.l^ uses a non-CGC framework where only the first power 
corrections to leading-twist collinear factorization are considered. 



5. Bulk features of multiparticle production in A+A, p+A and 
p+p collisions 

5.1. Multiplicity distributions in p+p and p+A collisions 

In recent years calculations of single-inclusive multiplicities within the CGC frame- 
work have been extended to multiplicity distributions and to their evolution with 
energy and system size. Much of this initial progress on multiparticle production 
in the non-linear regime is based on the framework of classical fields. Some basic 
insight into quantum evolution effects on various moments of the multiplicity dis- 
tributions has also been developed^^ but an explicit discussion of their properties 
is still lacking. 

Assuming dominance of classical fields, the probability to produce q particles is 
given 

(c'yi ■■■'^2/9) conn. """^ (c^yi) {dyg^} ^^^^ 
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C,^'^^ (15) 



This expression is valid with logarithmic accuracy and was derived under the as- 
sumption that all transverse momentum integrals over pT,i ■ ■ -PT.q are effectively 
cut off in the infrared at the non-linear scale. 

The fluctuation parameter k in eq. ([TSl) is of order 



k^{N!~l)QlS^ . (16) 

The numerical prefactor (in the classical approximation) has been determined by a 
numerical computation to all orders in the valence charge density in ref.^'^' 

The multiplicity distribution is therefore approximately a negative binomial dis- 
tribution (NBD),§2j 

_ T{k + n) n^k^ 
^'^^ ^ V{k)T{n+l) (n-Hfc)«+fc ■ ^ ' 

Indeed, multiplicity distributions observed in high-energy p-l-p collisions (in the 
central region) can be described reasonably well by a NBD, see for example refsl^^llsll 
The parameter k^^ determines the variance of the distributioiHI and can be obtained 
from the (inclusive) double-gluon multiplicity: 

rf2/idy2/co„„. k \dyi/ \dy2/ 

This expression shows that the perturbative expansion of k^^ starts at 0{a'^) since 
the connected diagrams on the l.h.s. of eq. (IT51) involve the same number of sources 
and vertices as the disconnected diagrams on the r.h.s. of that equation.l^ Since 
the mean multiplicity in the classical limit is of order n = (dN/dy) ~ 1/as it 
follows that fi/k ^ l/ctg ^ 1- This, in fact, corresponds to the limit where the 
NBD multiplicity distributions exhibit KNO scaling.!^ The semi-classical strong- 
field limit of a Gaussian effective theory relates the emergence of KNO scaling in 
p_L-integrated multiplicity distributions to properties of small- a: gluons around the 
saturation scale .1221 

KNO scaling conjectures that the particle multiplicity distribution in high- 
energy hadronic collisions is universal (i.e., energy independent) if expressed in 
terms of the fractional multiplicity z = n/h. This is satisfied in first approximation 
in the central (pseudo-) rapidity region at energies > 200 GeV as shown in fig. [S] 
For a detailed review of pre-LHC data see ref.l^ 

However, more precise NBD fits to the multiplicity distributions reveal that n/k 
increases somewhat with energy, and so KNO scaling is not exact. An increase of 
n/k ^ 1/as with energy may be partially explained by running of the coupling 

"^The width is given by n \/k^^ -\- n^^ ~ fij yjk where the latter approximation applies in the limit 
fijk > 1. 
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Fig. 5. KNO scaling plot of charged particle multiplicity distributions in non-single difTractive 
PP / PP collisions at various energies as measured by the UA5p3 ALIClP^ and CM!^^ collabora- 
tions, respectively. Note that we restrict to the central rapidity region I?;! < 0.5 and to the bulk of 
the distributions up to 3.5 times the mean multiplicity. 



with Qs{x) but this goes beyond the classical approximation. Corrections beyond 
the classical limit are also indicated by the fact that at lower energies higher-order 
factorial moments Gq of the distribution are significantly different from the reduced 
moments CgP^ 



a^^inin-l)...in-,+ l)) M. (19) 

Note that the difference of Gq and Cq is subleading in the density of valence charges 
p. In fact, while eq. P3|) has been derived for a quadratic MV-model action which 
applies near infinite density, at finite density the action contains additional terms 
~ p4j23l[99| These operators provide corrections to the moments of the mutliplicity 
distribution.'22l 

Finally we also point out that approximate KNO scaling has been predicted to 
persist even for min-bias p+Pb collisions at LHC energies,-^ with deviations from 
the KNO curve from p+p collisions at a level of "only" ^ 25%. This is in spite 
of Glauber fluctuations of the number of target participants which of course exist 
for a heavy-ion target. The multiplicity distribution in p-l-A collisions thus tests 
the expectation that intrinsic particle production fluctuations are dominated by the 
dilute source, i.e. that k ^ min(Tyi, Tg) is proportional to the thickness of the more 
dilute collision partner. Furthermore, it can provide valuable constraints for models 
of initial-state fluctuations in A+A collisions which are currently of great interest 
(see section [5^ . 



5.2. The initial state in A+A collisions 

The CGC formalism also offers a framework which allows one to determine the 
initial state of a heavy-ion collision. This refers to the distribution of produced 
particles in momentum as well as transverse coordinate space which determines 
the initial condition for the subsequent thermalization stage, followed by (viscous) 
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hydrodynamics. Understanding the thermahzation process in high-energy heavy- 
ion coUisions is a topic which is currently under intense investigation; we refer the 
interested reader to a recent review by F. Gehs published in this journal!^ 

In practice, the "detailed" dynamics of thermahzation is often skipped and the 
pr-integrated distribution of produced particles in the transverse plane is used di- 
rectly to initialize hydrodynamic simulations. This appears reasonable if one is 
interested mainly in the features of dN/dy(PrT over length scales larger than the 
thermahzation time rth- On the other hand, density fluctuations on shorter scales 
could be washed out, see e.g. ref.'sSl for a concrete example applying classical Yang- 
Mills field dynamics. 

Instead of resorting to the CGC formalism a much simpler shortcut to the ini- 
tial dN /dydPrr is to assume that the local density of particles produced at a point 
tt is simply proportional to the average density of projectile and target partic- 
ipants, dN/ dyd^VT ^ {Pparti^T) + Ppart(''r))/2. This "wounded nucleon" model, 
often times also labeled "Glauber model" , clearly must capture the rough features 
of dN/dyd^TT- On the other hand, it is evident from the increase of (l/A'part) dN/dij 
or (l/A'part) dET /drj from peripheral to central Pb-|-Pb collisions by about a factor 
of two (see figE]) that such a model is far from accurate. There is no reason to be- 
lieve that the second harmonic moment of the density distribution in the transverse 
plane, the so-called eccentricity e = {y'^—x'^)/{y^+x^) would do better (the same 
goes for higher moments). Indeed, a variety of CGC models which do describe the 
centrality dependence of dN/dy predict higher eccentricity e than the "wounded 
nucleon model" .^^"SSHim] 

The "wounded nucleon model" for soft particle production can be improved 
by adding a semi-hard component. It has to incorporate impact parameter and 

dependent shadowing in order to smoothly interpolate to p-l-p collisions in very 
peripheral collisions and to ^ A^coU scaling at high p^j^ also, the energy dependence 
of the low-pT cutoff required by leading-twist calculations needs to be fixed carefully 
to reproduce measured multiplicities. 

The initial spatial particle distribution exhibits large fluctuations. They manifest 
in non-zero elliptic flow V2 in central heavy-ion collision4i221 as well as in a large 
"triangular flow" component ^ig.liMlQne source of fluctuations is due to the locations 
of participant nucleonsfi23J these have also been incorporated early on in Monte- 
Carlo implementations of the fcr-factorization formula with KLN UGDs.'^'^ 

However, even for a fixed (local) number of participants there are intrinsic par- 
ticle production fluctuations. This is most evident from the wide multiplicity distri- 
bution in non-single diffractive p-fp collisions (see section [Q]) . A suitable extrapo- 
lation to A+A collisions has to be included both in "wounded nucleon"'i2^ as well 
as in CGC based^ Monte-Carlo models. In the CGC approach, intrinsic particle 
production fluctuations are expected to occur on sub-nucleon distance scales on the 
order of ^ l/QsW^&s shown in flg.[6l It is interesting to note that early CGC initial 
state models which do not incorporate intrinsic particle production fluctuationd^HIini 
appear to be inconsistent with the distributions of angular flow harmonics measured 
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by the ATLAS coUaboratiorP^ while recent approaches can describe higher eccen- 
tricity harmonics rather wehP^ The "cosmology of heavy-ion collisions" will be to 
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Fig. 6. Left: fluctuating energy density distribution at midrapidity and time t = 0.2 fm/c in an 
A-j-A collision.13^ Right: Scaled distributions of various azimuthal harmonics for Pb+Pb coUisions 
at 2.76 TeV (20%-25% centrahty class) from the CGC implementation of ref.C22l compared to 
ATLAS data. 



understand the nature, scale, magnitude and evolution of these fundamental QCD 
fluctuations in detail. It represents a very exciting avenue for future research. 

5.3. Multiplicities and transverse energy 

We now turn to the centrality dependence of the multiplicity and transverse en- 
ergy in heavy-ion collisions at LHC energies as obtained from A:T-factorization with 
rcBK UGDs. We recall that fcT-factorization is not expected to provide accurate re- 
sults for pt integrated observables in A-f A collisions .'^'^ Furthermore, we restrict 
to the initial "gluon liberation" process and do not account for subsequent gluon 
multiplication towards thermalization.^'^ ^'^^ Nevertheless, it is of course impor- 
tant to check that basic trends such as the centrality dependence of dN/dy are 
consistent with theoretical expectations regarding the dependence of the saturation 
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momentum on the thickness of a nucleus P^l^^ 
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Fig. 7. Left: Ccntrality dependence of the charged particle multiplicity at midrapidity, rj = 0; 
Pb+Pb collisions at 2760 GeV. We compare our calculation for two UGDs to data by the ALICE 
collaboration. Right; Centrality dependence of the transverse energy at r? = 0. 



We shall focus on the centrality dependence of the charged particle multiplicity 
at central rapidity, 77 = 0, which we determine along the lines described at the 
beginning of section S) The result is shown in Fig. [7] for the UGDs with MV-model 
(7 = 1) and AAMQS (7 = 1.101 and Ql{xo) Ta or Ql{xo) T^'') initial 
conditions. We use K = 1.43 for the former and K — 2.0, 2.3 for the latter. The 
number of final hadrons per gluon is Kg = 5 in all three cases. 

Aside from normalization factors, all UGDs give a rather similar centrality de- 
pendence of the multiplicity, and are in good agreement with ALICE data.'^^^'^ii' 
On the other hand, they differ somewhat in their prediction for the transverse en- 
ergy. This is of course due to the fact that the 7 > 1 initial condition suppresses 
the high-fcf tail of the UGeE|. With the iiT-factors mentioned above these UGDs 
match the measured Et in peripheral collisions. This is a sensible result since one 
does not expect large final-state effects in very peripheral collisions. For the most 
central collisions the energy deposited initially at central rapidity is about 0.5% of 
the energy of the beams and exceeds the preliminary measurements by ALICEfii^ 
and CM^SUby roughly 50%. This leaves room for —p AV work performed by lon- 
gitudinal hydrodynamic expansion.lii^ltii^ 

5.4. The CGC and its ridge 

In this section we discuss briefly the "ridge" structure observed in heavy-ion col- 
lisions at RHIC and in high- multiplicity p+p and p+Pb collisions at the LHC. A 
review dedicated entirely to this phenomenon has appeared recently in ref.'^^ 

"^One should keep in mind though that our estimate of the initial transverse energy carries a 
significant uncertainty of at least ±15% related to our choice of if-factor; it is not determined 
accurately by the multiplicity since the latter involves only the product of /C-factor and gluon — > 
hadron multiplication factor Kg. 
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The "ridge" refers to two-particle correlations at sma/Z relative azimuth, A(f> ^ 0, 
which extends over at least several units of relative (pseudo-)rapidity A77, perhaps 
even all the way to the sources. These do not correspond to particles produced in 
a binary collision since those would appear in opposite azimuthal hemispheres, due 
to transverse momentum conservation. 

What makes the ridge particularly interesting is that particles separated by 
a large pseudorapidity interval A77 are causally disconnected and hence can not 
be correlated, unless they were produced early on at a time on the order of 
T ^ Tf exp(— A77/2) or less.^^^ This argument relies on the assumption of Bjorken 
longitudinal boost invariance. 

The ridge has been first observed in Au+Au collisions at RHIC.'^i^'i^ The 
"heavy-ion ridge" is difficult to interpret from a CGC perspective since both the 
amplitude as well as the azimuthal coUimation of the correlation are certainly mod- 
ified by final-state interactions with the dense medium; those might be interpreted 
as jet-medium interactions, hydrodynamic flow etc. 

The discovery by CMS of similar ridge correlations in high multiplicity p-l-p 
collisions at 7 TeV therefore represented a major breakthrough. The duration of any 
final-state interactions for such small systems (in the transverse direction) clearly 
must be much shorter than in heavy-ion collisions: hydrodynamics predicts that 
(energy) density gradients decay on a time scale ~ R/cs where Cs — I/a/S denotes 
the rclativistic speed of souncl^. On the other hand, within the CGC framework it 
is only the local density of gluons that matters but not the size and life-time of the 
system; the ridge correlation is explained as an initial-state effect, anyhow. 

CMS has recently confirmed such ridge correlations also for high-multiplicity 
p-f Pb collisions at 5 TeV as shown in fig. [H] While the radial extent of the collision 
zone and the particle density is comparable to that from high-multiplicity p-l-p 
collisions the observed amplitude of the ridge is much higher. One may anticipate 
that this should be related to the fact that multiplicity fluctuations in p+Pb are 
helped by the presence of Glauber fluctuations of A^part we shall return to this 
point below. 

The origin of the correlation can then the explained in two ways. In the Kovner- 
Lublinsky pictur^SllESlEIIl ^^^^ incoming gluons from the projectile wave function 
with different rapidities scatter off the dense target field. If their relative impact 
parameter is less than the correlation length 1/Qs,t of the dense gluons in the target, 
and if their color charges are alike, they prefer to scatter to a relative azimuthal 
angle Acf) ^ and lead to a "near side" ridge. Projectile gluons with opposite color 
charges, a configuration which is equally likely to occur, on the other hand prefer to 
scatter to A ^ tt. Ref.'^^ makes the interesting prediction that this "degeneracy" of 
the near-side and away-side ridges is lifted in case of projectile quarks; this could in 

'Here, R denotes the scale of gradients which may even be substantially smaller than the overall 
radial extent of the collision zone. However, the radial expansion of very small density spikes is 
damped by viscosity. 
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Fig. 8. Left: the ridge in high-multiplicity p+Pb collisions at 5 TeV as seen by CMS p3 particles 
with 1 Ge V < pj' < 3 GeV. Right; ALICE ridge with subtracted jet yield reveal the dipole 
correlationP^SI ATLAS results for the p+Pb ridge are presented in ref(i^ 

principle be tested if one could trigger on hadrons which are more likely to originate 
from quark fragmentation. 

1 cflN 




t yAcl> 

TT 

Fig. 9. Illustration of the diagrams that dominate at different A(/>.E^ 

Another underlying picture for correlated production is as follows (see illustra- 
tion in fig. [3]). Two hard giuons from the same ladder are nearly back to back, by 
conservation of transverse momentum. On the other hand, those produced in sepa- 
rate ladders can have any relative angle. For correlated production, however, the two 
ladders need to be connected through the sources after squaring the amplitude. This 
is achieved by flipping the sources of the two ladders in the amplitude and its c.c; 
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the resulting diagrams are enhanced if the density of sources is non-perturbatively 
high.^^^ In fact, this leads to a variety of contributions (inch HBT-like correlations) 
as discussed in more detail in refP^ It has been realized early on that the "ridge" 
correlations are actually not uniform in but that the relative angle peaks at 
A0 ~ and A(j} ~ n^SMM 




Fig. 10. Left: correlated particle yield after ZYAM subtraction as a function of relative azimuth 
in different centrality bins. Right: associated near-side yield as a function of multiplicity N"^^'^'': 
dashed curves / symbols show how the yield increases with the number of participants in the Pb 
target, thin dotted lines track its dependence on the valence charge density in the proton. Both 
figures from refli^ 



The current state of the art quantitative results are due to Dusling and Venu- 
gopalan.E^Hmil They describe not only the small-angle but also the back to back 
correlations and in fact the entire "matrix" of correlations as a function of the 
transverse momenta of trigger and associated particle, and for different bins of mul- 
tiplicity (the latter shown fig. [TUl left). We do not delve into details here but do 
return to a point mentioned above. Fig. I10| right, shows how the associated yield 
increases with the overall multiplicity iV^°^'"°. For p+p collisions a very high multi- 
plicity is due to rare extreme density fluctuations in either of the nucleons; in p+Pb 
collisions high multiplicity fluctuations are helped by the presence of additional 
Glauber fluctuations of the target thickness. Folding the two, ref.E^H obtains a very 
natural explanation for the increased near-side yield in p+Pb as compared to p+p 
at the same muhiphcity A^°rk'"''- 



6. Conclusions 

Fourty years after the discovery of QCD, it is an established fact today that the 
density of gluons in a hadron grows rapidly with rapidity. Present day colliders have 
the potential to make us understand at what light-cone momenta and transverse 
distance or momentum scales non-linear color field interactions set in. Collisions 
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with large heavy-ion targets are the most promising environment for this physics. We 
hope we could shed some light in this review on the prospects of present p+A, A+A 
and possible future e+A programs to answer some related fundamental questions 
in QCD. The so-called "Color Glass Condensate" formalism to high-density QCD 
provides a systematic framework to address many of the underlying phenomena 
in a consistent way. Nevertheless, at present phenomenological applications often 
differ in terms of implementation and involve various approximations. Quantitative 
comparisons to data are therefore crucial in order to test the accuracy of our current 
tools such as BK evolution in the running-coupling approximation. The forthcoming 
p+Pb run at the LHC will advance tremendously our understanding of high-energy 
interactions and of the initial state of heavy-ion collisions. 
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